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We address the problem of how simple a solution can be for a given quantum local consistency instance. 
More specifically, we investigate how small the rank of the global density operator can be if the local constraints 
are known to be compatible. We prove that any compatible local density operators can be satisfied by a low rank 
global density operator Then we study both fermionic and bosonic versions of the A''-representability problem 
, as applications. After applying the channel-state duality, we prove that any compatible local channels can be 

y—^ • obtained through a global quantum channel with small Kraus rank. 

o: 

^ , PACS numbers: 03.65.Ud, 03.67.Mn, 89.70.Cf 

["T I . I. INTRODUCTION 

Understanding the various correlations and relationships amongst the different parts of a many body quantum system is one 
of the most difficult challenges in quantum theory. It is well known that the reduced density operators defined by partial traces 
characterize subsystems. Consider a system of three parties A, B and C: If all three two-particle density operators p"*^, p^'^ 
and p^'-^ are consistent with some global density operator p^^'-^, they must satisfy Ti-a{p^^) — Trc(p^'^), Ttb{p^^) ~ 
Trcip'^'^) and TtaIp'^'^) = TrB(/9^'-^). This is a necessary but not sufficient condition for the existence of p^^'^ . As the 
• particle number N becomes very large, the coiTelations between local subsystems become much more complicated. In general, 
^ ] local consistency is the problem of deciding whether a given collection of subsystem descriptions is consistent with some state 
Qi of the global system, or the problem of finding necessary and sufficient condition for consistency of subsystem descriptions. It 
is also called the quantum marginal problem in literature [1-3]. The community observed the relation between the spectrum of 
. bipartite quantum states and certain representations of the symmetric groups very recently. The consistency conditions for some 
^ ■ special classes of quantum states were then given in [1-4]. For general states with overlapping margins, the situation remains 
] unclear. 

■ If the particles under consideration are fermions instead of qubits, the local consistency problem has been known as the N- 
representability problem, which arose initially in the 1960's in connection with calculating the ground-state energies of general 
. interacting electrons [_ ]. 

" It was only recently shown in [6, 7] that both deciding the local consistency and deciding the local consistency for fermions 
, are QMA-complete, meaning both the consistency problem and the A^-representability problem are computationally at least as 
' hard as any other problem in the complexity class QMA. Here, QMA is the quantum analogue of the complexity class NP. 
I , Consequently, it is unlikely to have efficient algorithms for local consistency problems, even on a quantum computer. And very 
I> ■ recently, Wei, et al. proved that the bosonic version of the iV-representability problem is also QMA-complete[: ]. 

Though the local consistency problem is theoretically hard in the worst case, it is still worth exploring the potential solutions. 
There are various approaches scattered through the literature on this subject. Linden, et al. proved that almost every three-qubit 
] pure state can be uniquely determined among all states by their two-party reduced states['']. A related fermionic version was 
discussed in [10], whose results indicate that almost any three fermion pure state is uniquely determined among all states by 
their two-particle reduced states, though it was not stated explicitly. Linden et al.'s result was later generalized to A^-particle 
sy stems [1 1]. 

In this paper, we will focus on another direction of the local consistency problem. We are interested in how simple the 
solution can be, or more specifically, how small the rank of a solution can be. The same question for bipartite quantum system 
without overlapping margins was discussed in [ i ], but their approach seems technically difficult to generalize to the case with 
overlapping margins. In this work, we consider this problem regarding the rank of the solution in a very general setting, for 
multipartite quantum systems with overlapping margins. We provide a rank reduction based approach. We also show that some 
useful results from convex analysis can be applied to this problem directly, though it leads to a slightly weaker bound than 
ours[17]. Then we will apply our results to ferminonic and bosonic systems. Finally local consistency problem for quantum 
channels will be addressed. 

We now state our main result. For a given finite-dimensional Hilbert space H, B{H) will denote the space of bounded linear 
operators acting on H. For the n-fold tensor product 'H'^" and any integer i < n, Xi and Zi are general Pauli X-gate and general 
Pauli Z-gate for i-th qudit respectively. Formally, the local consistency problem can be stated as follows. 
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Problem 1. Consider a multipartite quantum system A1A2 ■ ■ - An with Hilbert space 'HaiA2-- An = T~Lai ®'H-A2 ®' ' '®'^An- 
Given a set of reduced density operators pi{i = 1, 2, • • • ,k) where each pi acts on a subsystem li of A1A2 ■ ■ ■ An, or li C 
{yli, yl2, • • • , An}- The local consistency problem is to address the existence of a global density operator p G BiJiAiA^-An) 
satisfying Tr/c p = pi for any 1 < i < k. If here is the complement of I i. 

An instance of the local consistency problem is a collection of pairs: a reduced density operator pi and corresponding 
subsystem li. For any given instance {(pi,/i)}i=i of the local consistency problem, if there is a global density operator 
p G B{7iAiA2 - A„) satisfying Tr/c p = Pi for any \ <i <k, then we say {(pi, Ii)}^=i is a compatible instance. 

In this paper, we will show that compatible local density operators are reduced states of some simple (low rank) global density 
operator. More specifically, we have following theorem. 

Theorem 1. For any compatible instance of local consistency problem {(pi, there is a solution with rank no more than 

(rankpi)^- 

We are primarily interested in instances where the whole system is "H®" and only no more than c-particle reduced states are 
known. In this case, the number of reduced states should be no more than (") , the rank is bounded by polynomial in n while the 
rank of a general density operator should be exponential in n. 

The paper is structured as follows. After introducing the requisite background material in section I, we give a proof of above 
theorem in section II. We then apply these ideas to the iV-representability problem in section III. In section IV we introduce the 
notion of consistency of quantum channels. Some examples are illustrated in section V. 

II. PROOF OF MAIN THEOREM 

Proof of Theorem 1. Since the local density operators of this instance are known to be compatible, we can start from some 
solution p, which is a density operator acting on 7iAiA2 - A„ and Tr/c p ~ pi for any 1 < i < k. Using a spectral decomposition 
of p we may write 

r 

where r = rank(p) and \tpi), \ip2), • • • , |'0r) are mutually orthogonal unit vectors. We will show that when r is large, then we 
can find another solution p' to the same instance and rank(p') < rank(p). 
Consider the spectral decompositions again and write 



P. = Pf IV'f > (V'f I for all 1 < z < A- (2) 

where r; = rankp^ and |V'2''')i ' ' ' ; IV'n'') are mutually orthogonal unit vectors for any i. 

Consider the following set: 

M = {X G B(supp p) : n,(Tr/o X)\l, = O for any 1 < i < k}. (3) 

Here, for each i, 11; is a projector on supp p^. 

Now, is a subspace of dimension at least — rank^IIi. 

If — X]i=i rank^IIi > 1, then is not empty. Let's say there is a non-zero X G M which implies X^ G A4 too. Thus 
both Hi = X + X^ and H2 = i{X — X^) are traceless Hermitian operators in M. Note that Hi and H2 cannot be zero 
simultaneously when X is not the zero operator. Without loss of generality, let us assume H = i?i(or H2) is non-zero. 

H is chosen from i3(supp p), or equivalently, there is some e such that p ± eH > which follows pi ± e Trjc^H) > 
for any i. Thus Hermitian operator Tti<:{H) lies completely in i?(supp pi) that implies ni(Ti7<: H)Ili = Trjc H and then 
Tri? = TrTr/eif = 0. 

Since the operator H 0,Tt H — contains both positive and negative eigenvalues, the same holds for p — XH for A >> 1. 
Hence there exists an intermediate value < A < 00 for which the operator p — XH is nonnegative, but not strictly positive, i.e. 
p — XH is a degenerate density matrix we are looking for. 

Now, for any solution p to a given instance {pijli}^^^, if — X]!Li rank^IIi > 1, we can always find a non-zero traceless 
Hermitian operator H G A^,and then another solution p' = p — XH to the same instance with rank less than raiikp . Thus , 
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by repeating this procedure until above quadratic inequality doesn't hold anymore, we will finally end with a solution cr with 

ranker < [\JYh=i (rankpi)^]. □ 

Corollary 1. For any instance of the local consistency problem with given local systems if the solution set is nonempty, 

then there is a solution with rank no more than \_'\J^\^i (dim/i)^] . 

Remark 1. Barvinok proved that if there is a positive setnidefinite matrix X satisfying 

{Q„X)=q„ Vl<^<fc (4) 

where Qi, • • • , Qk are symmetric matrices and qi, - ■ ■ ,qk ore complex numbers, then there is a positive semidefinite matrix X* 
satisfying the same equation system and additionally rank X* < [^^-^^—^J [12 ]. The main ingredients of his proof are the duality 
for linear programming in the quadratic form space. After applying Barvinok' s theorem to the local consistency problem, we 

will have a similar rank reduction which will lead to a solution with rank no more than \J 2 X]i=i dim^/i. Thus this result is 
weaker than ours. 



III. APPLICATION: Af-REPRESENTABILITY PROBLEM 



In this section, we will study the A^-representability problem, which is a fermionic analogue of the local consistency problem. 
The bosonic version of iV-representability is also addressed later. 
We first restate the A^-representabiUty problem as follows. 

Problem 2. Given a system of N fermions where each particle has d energy levels, and a k-fermion state p of size (^) x (J^), 
determine whether there exists an N-fermion state cr such that Trfc+i_... .Ar(cr) = p. 

According to the Pauli exclusion principle, no two particles can occupy the same state, thus we can always assume d> N. 

The space of A^-fermion pure states is mathematically described as the A^-th antisymmetric tensor product of with dimen- 
sion (^), denoted as A^C^. It is the span of all iV-fold antisymmetric tensor products of vectors xi, X2, • • • , xn in which 
is defined as 

1 ^ 

Xi A X2 /\ ■ ■ ■ /\ Xn = y £pXp(i) <Sl Xp(2) ® ■ ■ ■ ® Xp(M)- (5) 

Here, P goes through all permutations of N indices and ep is the signature of P. So £p is 1, if the number of even-order cycles 
in P's cycle type is even, and —1 otherwise. 

Similarly, the space of A^-boson pure states with d energy levels corresponds to the A^-th symmetric tensor product of Cd with 
dimension (^^~^), denoted as V^C^. 

For more information about A^-th symmetric/antisymmetric tensor product, please refer to [13]. 

For the A^-representability problem, there is a similar rank reduction as follows. 

Theorem 2. Suppose we are given a system of N fermions where each particle has d energy levels and a k-fermion density 
operator p of size (^) x (^). Assume there exists an N-fermion state a such that Trfc+i.... ,7v(cr) — P- Then there also exists an 
N-fermion density operator a' with Trfe+i.... ,Ar(cr') = p and ranker' < rank p < (^). 

The proof is similar to the proof provided in Section II, with minor modifications. Observe that the whole rank reduction 
in our approach is processed in supp(p). After introducing additional symmetry to the global system, the rank reduction also 
works by replacing ®^Crf with A^Cd or V^C^. 

Similarly, we will also get the following theorem for the bosonic version. 

Theorem 3. Suppose we have a system of N bosons where each particle has d energy levels and a k-boson density operator 
p of size if'^^r^") X {^^^^\ Assume there exists an N-boson state a such that Trfc+i^... ,jv(cr) = A then there also exists an 
N-boson density operator a' with Trfe+i.... ,Ar(cr') — p and ranker' < rank p < ('^~''^~^) • 



IV. LOCAL CONSISTENCY PROBLEM FOR QUANTUM CHANNELS 



In this section, we will investigate a new type of consistency - the consistency of quantum channels. A quantum channel 
is a device which transmits classical bits or quantum states. Mathematically, it is a linear map which maps any quantum state 
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on some Hilbert space 'Hi to another state on some Hilbert space 7^2- Furthermore, a quantum channel can be described by a 
completely-positive, trace-preserving map <I>. 

Generally, for a quantum channel from system Hi to H2, we shall think of Hi as part of a closed composite system Hi (E) H'l 
and H2 as part of another closed composite system H2 ® H2 which has same dimension as Hi <E) H'l- Therefore, the evolution 
from Hi (E) H'l to H2 E) H'2 can be described by some unitary operator U . The quantum channel $ is then described as 

$(p)=Tr«.([/(p®|0)(0|«;)[/t). (6) 

By Stinespring's dilation theorem on completely positive maps, $ must take following form 

N 

4=1 

where K[s are some operators, called Kraus operators of $. Trace preservation of $ is equivalent to the sum X^iLi ^l^i 
equaling the identity operator The number of Kraus operators N is no more than dim "Hi dimH2, and the minimum number of 
N is called the Kraus rank of $. In some sense, the smaller the Kraus rank is, the simpler the channel is. 

The concept of channel consistency is quite intuitive. Consider the following scenario, in which there is a channel from some 
large system Hi to some large system H2- Here, the mapping is from B{Hi) to B{H2)- A local observer Alice can only gather 
information from part of Hi, say Hi and part of H2, say H2', therefore, she has information about the partial mapping from 
B{Hi) to B{H2)- Another observer Bob has his information about the partial mapping from some B{Hf) to some _B('H^). 
So do any other observers. Several questions naturally arise. How much information about the global quantum channel can 
be known when Alice, Bob and other observers disclose their local information? If every observer has a description of some 
partial mapping, is there a global channel satisfying all these local constraints? Can we find some simple channel satisfying local 
constraints if the local descriptions are known to be compatible? 

There are some subtle differences between state consistency and channel consistency. The most confusing part is, how to 
describe part of a quantum channel, or the partial mapping from a local system to another local system? There is no doubt that 
part of a quantum state p, described by applying some partial trace on p, is definitely positive-semidefinite and trace 1. Therefore, 
part of a quantum state is again a quantum state. However, in the quantum channel setting, the analogue of the above property 
is not so straightforward. One may even doubt part of a quantum channel may not be a quantum channel at all. Thanks to the 
channel-state duality, we can define sub-channel of some quantum channel as the following. 

Given any channel ^' : B{Hi) — >■ B{H2), we can always write a corresponding state of '5 to be 



""^dh^i ^ (8) 

p,q=l 

where is an orthonormal basis of Hi- The state axj, is called the Choi-Jamiolkowski state of 4* and the association 

above defines an isomorphism between linear maps from B{Hi) to B{H2) and operators in B{Hi ® H2), called the Choi- 
Jamiolkowski isomorphism. Its rank is equal to the Kraus rank of cr*. 

Therefore, for any quantum channel : B{Ha <E)Hb) ^ B{Ha' ® Hb ), we can define channel : B{Ha) ^ B{H\) 
by taking the reduced density operator of the Choi-Jamiolkowski state tr^ as its Cho-Jamiolkowski state a^^A . 

Observe that 

^ Aim'H.AB 

= Tva{{^—^ TrB|p)(g|®TrB'Vl'(|p)(q|))(p®/^0) (9) 

p,q=l 

diiTiT-L AB 

= ^--^ E ^rAB{{\p){q\{p(ElB))(ETTB'^{\p){q\)) (10) 

= TrB^ip®-^-), (11) 
dim Hb 

"ii^ acts exactly the same as does between B{Ha) and B{H'y^). Hence we call sub-channel of \E' from B{Ha) to B{H'j^). 

By adopting above definition, we will address the following question: how simple the global channel can be, or more specif- 
ically, how small its Kraus rank can be if the sub channels are known to be compatible. Mathematically, the local consistency 
problem for quantum channels can be stated as follows. 

Problem 3. Consider two multipartite quantum systems A1A2 ■ ■ ■ An and B1B2 ■ ■ ■ Bm with Hilbert spaces HaiA2 - A„ = 
Hai ® T~La2 ® ■ ■ ■ ® Ha^ '^^d T~LBiB2-Bm — ^Si ® ® ■ ■ ■ ® T~LBm respectively. Assume a set of local quantum channels 
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{$; : Z = 1, 2, • • • , fc} ;s given. Each $; maps states on Hilbert space H/, = ®Ai£ii T~LAi to states on Hilbert space "Hj^ = 
®Bj^Ji where Ii C {^i, ■ ■ ■ , An} and Ji C {Si, B2, • • • , Bm}- The local consistency problem for quantum channels 
is to address the existence of a quantum channel : B{'HaiA2-A„) ^ B{Hq-^j3^...q^) satisfying ^i{p) = Tr jc ^{p(S> ^[jj-^^^ ) 
for any p G D{'Hii) for any 1 < ' < fc- 

By taking the Choi-Jamiolkowski states of each Problem 3 can be converted to the existence of the global density operator 
Then we can apply Theorem 1 to the quantum system A1A2 ■ ■ ■ AnBiB2 ■ ■ ■ Bm with Hilbert space 'HAiA2 - Ar,BiB2---Bm = 
Hai €5 'Ha2 ® • • • $5 HAr, ^ ^Bi <8) 7^62 (8) • • • (8) Hsm- Since a set of reduced density operators cr$, (/ = 1, 2, • • • , fc) is given 
and each acts on a subsystem /; J; C {Ai,A2, ■ ■ ■ , An, Bi,B2, • • • , Bm}. According to Theorem 1, we can always find a 

density operator a satisfying the local consistencies and rank a < yjj^i (dim/; dim J; )^. 

Theorem 4. Assume there exists a global quantum channel $ : B{'HaiA2 - A„) ~> 5('HeiB2 •■B,„) satisfying = 

Trjc $(p (g) j.^^ ^ ) for any p G D[Ti.i^) for any 1 < I < k. Then there also exists a quantum channel $' satisfying the 
' 'i 

same local constraints such that $' can be expressed with no more than y X); (dim/; dim Jif' Kraus operators. 



V. SOME EXAMPLES 

Example 1. Consider an n-qubit quantum system A1A2 ■ ■ ■ An with Hilbert space 'HaiA2 - A„ = Cf"". We are interested in 
the n-qubit states p such that any k-qubit local density operator of p is 
Obviously, p — ^ is a trivial candidate with the maximal rank 2". 

From Theorem 1, there exists some n-qubit density operator p satisfying same local consistency and rank(/5) € 0(2*^ (^))- 

As a corollary, when k = 2 is fixed, then rank(p) € 0{n). 

Indeed, for k = 2, there are always pure state (i.e. rank ~ 1} solutions for any n > b. One such example could be a graph 
state [14, 15] Itjjn) on a ring, which is a common eigenstate of eigenvalue 1 of the Pauli operators gi = {Zi^iXiZi^i} for 
i = 1, 2, . . . , n, where Zq = Zn, Zn+\ = Z\. That is, |V'n) = IV'n) for i — 1, 2, . . . , n. Note that 



Pn = \i^n){^n\ = ^J[{I + 9^)■ (12) 

It is then straightforward to see that any 2-local density operator of the n-particle state pn is 

In general, for any fixed k, there do exist n-qubit graph states such that any k-local density operator of the graph state is ||, 
for large enough n [16]. 

Example 2. Consider a system of N bosons where each particle has 2 energy levels. The 2-boson maximally mixed state is 
defined as 

= i(|00)(00| + |11)(11| + (13) 
Obviously there exists a non-degenerate N-boson maximally mixed operator M^j^^ such that Tra,... ^n{AI^^) ~ ^J^b ^ ■ Then it 

(2) 

follows from Theorem 3 that there exists another N-boson density operator a with Trs^... ,Ar(cr) = and ranker < 3. 

Wfe can choose < P < ^^T^ '^"'^ '^v 

Notice that 

|«lJ2---«Ar) (15) 

il+i2H ViN=V 

= ioo> \h---iN) + m) + \m Yl h---iN) + \n) y h---iN), m 

isH hiN=P isH \-iN=p-l isH hijv=p-2 
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and then 



Trs,... ,Ar(crp) 



(17) 





01 + 10, ,01 + 10 



|01 + 10)(01 + 10| + ( 



2jV-3p+l (pj) 
6(7V-p) 



(18) 



^/2 ^/2 



I). 



(19) 



Therefore, {cp} n-i <p< 2n+i ;s a family of N -boson density operators with rank 3 a«c/ every 2-local density operator of any 
CTp ii f/ie bosonic maximally mixed state . Furthermore, when N = l(mod3), we will have rankfcr w-i ) = 2. 



In this paper we addressed the problem of how simple a solution can be for any given local consistency instance. More 
specifically, how small the rank of a global density operator can be if the local constraints are known to be compatible. We 
provided a reduction based approach to this problem and proved that any compatible local density operators can be satisfied with 
a global density operator with bounded rank. Then we studied both fermionic and bosonic versions of the A^-representability 
problem as applications. After applying the channel-state duality, we proved that any compatible local channels can be satisfied 
with a global quantum channel which can be expressed with a small number of Kraus operators. 

This paper represents a preliminary step toward understanding the structure of solutions to the local consistency problem. 
There are many open questions from this approach deserving further investigation. For example, though the local consistency 
is known to be QMA-complete in general, efficient algorithms are still possible for some classes of instances. Since we know 
now the existence of solutions is equivalent to the existence of simple solutions, we can ask if it is possible to find more efficient 
algorithms for these classes? Further, we could ask if only spectra or other descriptions are known for subsystems, how simple 
can a solution be? 
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